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Abstract 

The purpose of this paper is to describe asymptotic formulas for determinants of 
a sum of finite Toeplitz and Hankel matrices with singular generating functions. The 
formulas are similar to those of the analogous problem for finite Toeplitz matrices for 
a certain class of symbols. However, the appearance of the Hankel matrices changes 
the nature of the asymptotics in some instances depending on the location of the 
singularities. Several concrete examples are also described in the paper. 



1 Introduction 



In the theory of random matrices jnj one is led naturally to consider the asymptotics of 



determinants of Fredholm operators of the form I + W + H where W is a finite Wiener-Hopf 
operator and if is a finite Hankel operator. This problem arises when investigating the 
probability distribution function of a random variable thought of as a function of the eigen- 
values of a positive random Hermitian matrix. The random matrix connections show that 
the constant term in the asymptotic expansion of determinants is fundamentally connected 
to the mean and variance of the distribution function. We will not describe the random 
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matrix connections any further, but simply refer the reader to |16| for general information 



and also [l], [|, H for more specific tie-ins to the random variable problem. 

The focus of this paper is to study the discrete analogue of this problem. This is not 
exactly the desired situation for those interested in random matrix theory. However, it is a 
natural starting place for cases where the random variable is discontinuous, since then the 
discrete nature of the computations make things a bit more accessible and the mathematical 
questions that arise are quite interesting in themselves. 

The discrete analogue of the Fredholm determinant problem is precisely to find an asymp- 
totic expansion of the determinants of the following Toeplitz + Hankel matrices 

M n (0) = r n (0) + # n (0). (1) 

Here the n x n Toeplitz and Hankel matrices are defined as usual by 

Tn(4>) = (0i-fe)",fcio> H n(4>) = (<Pj+k+l) n ~ k l - (2) 

The entries <pk are given by the fc-th Fourier coefficient of where is a sufficiently well- 
behaved function on the circle. If is continuous, even, and sufficiently smooth, then the 
continuous analogue of the problem (i.e. the Toeplitz + Hankel matrices are replaced by finite 
Wiener-Hopf + Hankel operators) is solved in [|J . There it is shown that the asymptotics are 
very similar to the ones given in the Szego-Kac-Widom Strong Limit Theorem. Indeed, it 
is only in the constant, or third order term that the answers differ. This is no surprise since 
if is continuous, then the Toeplitz operator is perturbed by a compact Hankel operator. 
However, if the symbol is singular, then the problem, as in the Toeplitz case, is not easy 
to solve. 

The purpose of this paper is to compute the asymptotics of detM n (0) as n — > oo for 
certain piecewise continuous functions 0. The main general result that we will obtain is as 
follows. We consider piecewise continuous functions of the form 

R 

0(e^) = 6( e <tf )t Af ( e *)t /3L ( e 4 ^))JJ*fr(e i( *- flr) ), (3) 

r=l 

where 

tp(e ie ) = exp(i/3(0-7r)), < < 2tt, (4) 

and b is a smooth nonvanishing function defined on the circle with winding number zero. We 
also need conditions on the parameters /?+, /?_, Pi, . . . , (3r. These conditions on the "size" of 
the jumps and the precise smoothness conditions on b will be described later on. In addition, 
we have to assume that 6\, . . . , Or G (—it, 0) U (0, n) are certain distinct numbers satisfying 
r + S 7^ for each r and s. The latter condition excludes piecewise continuous functions 
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with jumps at both a point on the unit circle and its complex conjugate. However, the 
function may have jumps at the points 1 and —1. Our conditions on b guarantee that the 
following functions 

M*) = exp ^t fc [log6] fc j , teT, (5) 

b-(t) = exp ^|>->g&]- fe j , teT, (6) 

are well defined and smooth. Here [log b]k denotes the k-th Fourier coefficient of the logarithm 
of b. Then the asymptotic formula reads 

detM„(0) ~ G[b] n n nM E M (7) 

as n — > oo, where 

G[b] = exp[log6] , (8) 

o 3^ ft. 3/3 2 /3_ A 

r=1 



£ M = E[b]F[b] 

x G(l + - P+)G(l/2 - (3 + )G(l/2)-\2n) l3+/2 2 3/3 + /2 

x G(l + - P-)G(3/2 - ^)G{?,/2Y\2nf-/ 2 2 3/3 -/ 2 



R 

ie sff/2+(3 r /2 , iB sff/2-f3 r /2 



x H G(l + (3 r )G(l - ft) (1 - e -«»r^r/2+ft-/2 ^ + e -i9 r j 

r=l 

x 6 + (l) 2/3 +6_(l)^+6 + (-l) 2/3 -6_(-l)^- 2 3/3+/3 



r=l 



X 
r=l 



JJ (1 _ e -^) 2 ^ (i _ e *r)P + * (! + e -w,)^-^ (! + e Wry 

r=l 

x n I 1 - e^"^)^' (! " e*^- 9 '))^' (l - e -(W)^ . (i ) 

l<s<r<_R 

The constants £"[6] and F[6] are defined by 

E[b\ = exp ^fjA;[logfe] fc [logfe]_ fc j , (11) 

F[b] = (^) 1/2 -p(-^f;Mlog^j, (12) 



and G(*) is the Barnes G-function p], |18| defined by 

oo 

G(l + z) = (2nY^e-^ z / 2 -^ z ^ 2 l[^l + z/k) k e- z+z2/2k ^ (13) 



k=l 



with being Euler's constant. 

It is interesting to compare this asymptotic formula with the corresponding formula for 
Toeplitz determinants. The asymptotic expansion of Toeplitz determinants for a certain 
class of singular generating functions is described by the Fisher-Hartwig conjecture. For 
instance, it is well known || H that if <fi is of the form 

R 

J8\ 



r=l 



b(e t9 )i[tp r (e^), (14) 



where b is a sufficiently smooth and nonvanishing on the unit circle, 9i, . . . ,9r G (— 7r,7r] 
are distinct numbers, and if |Re/5 r | < 1/2 holds for each 1 < r < R, then the asymptotic 
behavior of the Toeplitz determinants is given by 



detT n (0) ~ G[b} n n UT E T (15) 

as n — > oo, where 

n T = -f> 2 , (16) 

r=l 

R 

E T = E[b] ]jG(l+(3 r )G(l-(3 r ) 

r=l 

x Jj6 + (e^) ft '6_(e ie '')- /3 '- JJ - e*^'-^)^'. (17) 



A general account of the Fisher-Hartwig conjecture can be found in [10| and more recent 
work in |l^, It has been proved in many cases and reformulated in others. 

The paper is organized as follows. In Section we compute the asymptotics of det M„(0) 
in the case of smooth nonvanishing functions </> with winding number zero. As in the contin- 
uous analogue, the asymptotics are very similar to the Toeplitz case. In fact, we prove that 
the quotient detM n (0)/ detT n (0) converges to a nonzero constant. 

In Section |3| we recall several operator theoretic results, in particular those related to 
Toeplitz operators and Toeplitz + Hankel operators. 

In Section ^ we prove the asymptotic formula (]?]) in the special case of piecewise continu- 
ous functions (H) without jumps at 1 and —1 (i.e. with /3 + = /3_ = 0), without jumps at both 
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a point on the unit circle and its complex conjugate (i.e. 9 r + 9 S ^ 0) and under the assump- 
tion |Re/3 r | < 1/2 , 1 < r < R. As before we show that detM n (0)/ detT n (0) converges to 
a nonzero constant (althought the corresponding Hankel operator is not compact). In other 
words, if the symbol has jumps in the "proper" locations, then the asymptotic behavior is 
again like in the Toeplitz case. Note that the condition on the location of the jumps is 
extremely important in the Toeplitz + Hankel contrasted to the Toeplitz case. 

In Section || we show that the quotient det M n (0-?/>)/(det M n (0) detM n (ip) converges to 
a nonzero constant for certain functions <ft and ip. This result allows us to localize at cer- 
tain points on the unit circle, in particular at 1 or —1. However, it is not possible to 
localize at a point on the unit circle and its complex conjugate. Thus the localization 
result reduces the asymptotics for general piecewise continuous functions to those for func- 
tions t p (e ie ) and t p {e i{ - 6 '^) with a single pure jump at 1 or —1 and to those for functions 
^/3+( e ^ 6 ' _e '' ' > )^/3 r r( e * ( ' e ' +e '' ' > ) w ith two pure jumps. Note that it is exactly this last class of func- 
tions for which we are not able to determine the asymptotics in general. 

In Section |6] we then consider the case of functions tp(e ie ) and tp(e^ 9 ~^). In this case, 
the jump at 1 or —1 on the circle changes the nature of the asymptotics in the second 
order term. The computations are based on the fact that the corresponding matrices are 
Cauchy matrices times certain diagonal matrices. Finally, this result in conjunction with the 
localization and the results of Section |] gives the above mentioned main general result (|^). 

In the last section, we illustrate with additional concrete examples. The first class of 
examples is for piecewise continuous functions with two jumps either at ±1 or at ±1 These 
functions are special cases where the parameters describing the jumps are connected with 
each other in some way. The significance is that they show that the one jump results, for 
functions with jumps at i and — i cannot be pieced together to obtain the asymptotics for a 
symbol that has jumps at both these points. It should be pointed out that this does work 
for the points 1 and —1 and this is confirmed by the examples. Note that one special case 
of these examples is a piecewise constant even function with jumps at ±i. 

The second class of examples in the last section is for the even functions 

u a (e ie ) = (2-2cos#) a and u a (e i(9_ir) ) = (2 + 2cos#) a . (18) 

These function are singular or zero at 1 or —1, respectively, and are particularly interesting 
since they represent a more general class of examples of even functions. 

For random matrix theory even functions are most important. It would be helpful in the 
future to extend these results to that case and also to the continuous analogue of Wiener- 
Hopf + Hankel operators. But we believe the present paper is a good start and leave the 
other questions to some other time. 
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2 Operator theoretic results and Toeplitz + Hankel 
determinants in the case of smooth functions 



In this section, we compute the asymptotic behavior of determinants of Toeplitz + Hankel 
matrices M n (0) in the case of smooth nonvanishing functions defined on the unit circle with 
winding number zero. What we exactly mean by smoothness will be explained shortly. In 
the first part of this section, however, we will recall certain operator theoretic results that 
will be needed later on. 

We first introduce the following linear bounded operators acting on the Hilbert space 
£ 2 = -^2(2^+) of one-sided square- summable sequences. Given G L°°(T), define 

M(0) = T(<j>) + H(<j>) (19) 

where the Toeplitz and Hankel operators are given by the infinite matrices 

T{4>) = (&- fc )JLo, Htf) = (<Wi)S =0 - (20) 

Note that the Hardy spaces H°° and H°° consist of those functions <fi G L°°(T) for which the 
Fourier coefficients (f) n vanish for each n < or n > 0, respectively. We also write 

0(e* e ) = fte*), (21) 

and call <fi even if <p = 0. Finally, we introduce the projection P n acting on £ 2 by 

P n (/0, /!,•••) = (/o,/l,-..,/n-2,/n-l,0,0,...). (22) 

Note that T n {<j>) = P n T(0)P n , H n (<j>) = P n H{<f>)P n and M n (<j>) = P n M(<P)P n . 

It is well known that Toeplitz and Hankel operators are related to each other by the 
formulas 

= T(<f>)T(ij)+H(<f>)H($), (23) 
= T(4>)H(i>) + H(4>)T(i>). (24) 

If ip + G H°° and ip- G H°°, then we have 

T^_U + ) = T(^-)T(0)T(^ + ), (25) 
H^_U + ) = T(^)H(4>)T(i; + ). (26) 

Combining equations ( p3|) and (|24"D, it follows that 

M(#) = T(<j>)M(ip) + H{4))M$). (27) 
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This implies 

M{<j4) = M((j))M(ip)+H((j))M(^-ip). (28) 

If ijj is even, then the latter equation simplifies to 

M((j)ip) = M((p)M(ip). (29) 

These and other results concerning M(0) are discussed and proved in [0]. 

Two important notions are stability and strong convergence. Let A n be a sequence of 
operators. This sequence is said to be stable if there exists an n such that the operators 
A n are invertible for each n > n and sup n> HA" 1 ]] < oo. Moreover, we say that A n 
converges strongly on £2 to an operator A as n — > 00 if A n x — > Ax in the norm of £2 for 
each x G £2- When dealing with finite matrices A n , we identify the matrices and their 
inverses with operators acting on £2. It is interesting to note that stability is related to 
strong convergence of the inverses (and their adjoints) in the following sense. 

Lemma 2.1 Suppose that A n is a stable sequence such that A n — > A and A* n —>■ A* strongly. 
Then A is invertible, and A~ x — > A^ 1 and (A~ 1 )* — > (A^ 1 )* strongly. 

Another important set of operators is the ideal of trace class operators (see e.g. |15|]). For 
such operators, the trace "trA" and the operator determinant "det(J + v4)" are well defined 
and continuous with respect to A in the trace class norm. The following result shows the 
connection with strong convergence. 

Lemma 2.2 Let B be a trace class operator and suppose that A n and C n are sequences such 
that A n — > A and C* — > C* strongly. Then A n BC n — > ABC in the trace class norm. 

We proceed with describing the smoothness conditions. We therefore introduce certain 
function spaces. Let F£^ stand for the Banach space of all functions b G L X (T) for which 



\f4 




(i + M) 2/3 N 2 < 00, (30) 



1 /2 

and let W denote the Wiener algebra. It is well known that F£ 2 H W is a Banach algebra 
of continuous functions on the unit circle. The Besov class B\ is the Banach algebra of all 
functions b G L X (T) for which 

— / I b(e ix+iy ) + b(e ix ~ iy ) - 2b(e ix ) I dxdy < 00. (31) 

-7T V J-K 



Using results of Peller |T7| one can show that b G B\ if and only if both H (b) and H{b) are 



trace class operators. Moreover, the Riesz projection acts boundedly on B\. An equivalent 
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norm in B\ is given by |6q| + ll-^WIIi + ll^(^)l|i- Finally, one has the following continuous 
and dense embeddings 

Fl{ C B\ C (f^ 2 nW^j if > 1. (32) 

For more information on these and related classes of smooth functions we refer the reader 
to [|1(J and the literature cited there. 

The Besov class B\ exactly fits our purposes in the sense that the function h appearing in 
(0) will be assumed to be in B\. In order to simplify notation we denote by GqB\ the group 
of all nonvanishing functions in B\ with winding number zero. Remark that the asymptotic 
formula for Toeplitz determinants as given in ([15]) has been proved for functions of the form 
( fL4l) with b G GqB\ (see e.g. [0). The following proposition shows that all definitions 
involving the function b which were made in the introduction make sense. 

Proposition 2.3 Let b G GqB\. Then b possesses a logarithm log b G B\, and the constants 
G[b], E[b] and F[b] as well as the functions b + and b- make sense. Moreover, 

b(e i0 ) = b_{e i6 )G[b}b + {e ld ), < 6 < 2tt, 



is the (normalized) Wiener-Hopf factorization withb + G G(B\{~\H°°) andb- G G(B\C\H°°). 

Proof. Approximating the function b by polynomials and using that B\ is a Banach algebra, 
one can show that b possesses a logarithm log b G B\. Now one need only use the boundedness 
of the Riesz projection and the fact that B\ is contained in F£^ 2 D W. □ 
We now establish a limit relation for the quotient of two determinants with smooth 
generating functions. In Section [| this relation will be generalized to the case of certain 
piecewise continuous functions. 

Proposition 2.4 Let b G G B\. Then det M n (b)/ det T n (b) -> F(b) as n -> oo where F(b) 
is the operator determinant defined by 

F(b) := det (/ + T" 1 (0)#(0)) . (33) 

Proof. It is well known [0] that under the above assumptions on b, the sequence T n (b) is 
stable and the inverses converge strongly on £ 2 to the inverse of T(b). Since H(b) is trace 
class, we obtain that 

d ^(b) = dfltT » 1 ( 6 ) M »( 6 ) = det [p n + T-\b)P n H(b)P r , 

converges to the operator determinant F(b). Note that P* = P n — > / strongly. □ 
Next we establish the one main result of this section, the evaluation of the operator deter- 
minant F(b). The computation is especially remarkable since it relies on a "differentiation" 
trick, which was recently used in a modified form in M. 
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Theorem 2.5 Let b G G B\. Then F(b) = F[b]. 

Proof. Using the Wiener-Hopf factorization of b and formula ( |25| ) it is easily seen that 
T(b) = G[b]T(b-)T(b + ), and hence the inverse equals 

T-\b) = Glb^Tib-^TibZ 1 ). 

From equation (p6|) it follows that 

T-\b)H(b) = T(b+ 1 )T(bZ 1 )H(b-b+) = T- 1 (b+)H(b+). (34) 

Hence F(b) = F(b + ). Now let c = b + b + . One can conclude analogously that F(c) = F(b + ). 
Thus we are left with the evaluation of F(c). Obviously, [logc]„ = [log6]| n | for n G Z\{0}, and 
hence logc = logc. Since the Riesz projection is bounded on B\, we obtain that logc G B\. 
Next we define the Sj-valued analytic (in A) function 

c\ = exp(Alogc), AgC 

Note that the derivative of c\ with respect to A equals CAlogc. Let Y(X) be the analytic 
operator-valued function 

Y(X) = I + T-\c x )H(c x ) = T-\c x )M(c x ). 

It is easy to compute the inverse and the derivative of Y(X). We obtain 

F'(A)F- 1 (A) = -T- 1 (c A )T(c A log c)T- 1 (c A )M(c A )M~ 1 (c A )T(c A ) 
+ T-\c x )M{c x log c)M-\c x )T{c x ) 
= -T-\c x )T(c x \ogc)+T-\c x )M(\ogc)T(c x ). 



Note that c A = c A , and thus flS^) implies M(c A logc) = M(log c)M(c A ). Because det Y{X) ^ 
for all A, the scalar function y(X) = log det Y(X) is an entire analytic function. We conclude 

dX ~ detr " 1 

= tr (- T{c x \ogc)T~\c x ) + Af(logc)). 

Differentiating again yields 

= ^(-T^log^^-^c^+T^log^T-^c^T^log^T- 1 ^; 
= tr (-T(log 2 c) +T(logc)T(logc)) = - tr H (logc) H (logc). 
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The last equality follows from T~ l (c\) = T(c^\.)T(c^ 1 _) where c\ = c\-C\ t+ is the Wiener- 
Hopf factorization of c\. By repeated application of (|25f) , all occurring functions ca,± cancel 
each other. Hence the second derivative of y does not depend on A. Note that 

y(0) = logdet/ = 0, 

y'(0) = tr f-T(logc) +M(logcj) = trF(logc) 

since ca|a=o = 1. It follows that 

y(\) = -A 2 /2tr H (log c)H (log c) + A trif(logc). 
Because F(c) = detF(l) = expy(l), we obtain 

F(c) = exp ( - 1/2 tr H(log c)H (log c) + tr F(log c)) 

oo oo 

= exp ( - l/2^n[logc]„[logc]_ n + 1/2^ |[logc] n - (-l)"[logc] n }). 

n=l n=l 

Writing log c in terms of log b, we arrive at 

CO 

F(b) = exp(-l/2^n[logC + l/2 1og6 + (l)-l/2 1ogfo + (-l)). 

n=l 

This immediately implies the assertion. □ 
Finally, we can combine the previous results with the well known Szego-Widom Limit 
Theorem [[HJ. This limit theorem says that det T n (b) ~ G ! [6] n £'[6] as n — > oo for nonvanishing 

1/2 

functions b G -F£ 2 n W with winding number zero. 
Corollary 2.6 Let 6 G Go^ 1 . Then 

det M n (6) ~ G[6]"L;[6]F[6] as n -> oo. 

3 Further operator theoretic results 

The proofs of the results that will presented in the following two sections require further 
operator theoretic preliminaries. In particular, we need some results about Toeplitz operators 
and Hankel operators as well as about Toeplitz + Hankel operators M(<p) (see |TU[] and J7| 



for the general theory). First of all, in addition to the projection P n , we define Q n — I — P n 
and 

W n (f , /i, . . .) = (/ n _i, /„-2, • • • , A, /o, 0,0,...), 

V n {fa, fu ■ ■ ■) = (o,o,...,o,o,/ ,A,A,...), 

V-n(fo, fl, ■ ■ ■) = (fn, fn+1, fn+2, ■ ■ ■)■ 
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It is easily seen that W% = P n , W n = W n P n = P n W ni V n V- n = Q n and VL n V n = I. Note 
also that 

P n T(<j))V n = W n H(<P) } U_ n T(0)P n = H{<j>)W n . (35) 
Moreover, we have 

V_ ft #(0) = H(</>)V n , (36) 

WM^Wn = T n Q>). (37) 

Using equations ( pop we can write 

P n T(<p)Q n T(^)P n = P n T{(f)V n V^T^)P n = W n H0)H{,lj)W n . (38) 



Taking (^) into account we arrive at the fundamental identity due to Widom |19| 

T n {^) = T n (<f>)T n &) + P n H((j))H$)P n + W n H$)H(i;)W n . (39) 

The following result deals with strong convergence. For brevity of notation, we will 
henceforth write A = s — lim A n if both A n — > A and A* n —>■ A* strongly. 

Proposition 3.1 Let <p G L°°(T). Then 

T{4>) = s-lim T n (0), T(0) = s-lim W n T n (0)W n , 

M(0) = s-lim M n (0), T(0) = s - lim W n M n {<f>)W n . 

Proof. The relations without the W n 's are trivial because P* = P n I strongly. Also the 
second assertion is easy to show by taking account of (pTI). Finally, using fl35|) we obtain 

W n M n {4>)W n = T n Q>) + W n V- n T((f))P n = TM+P n T(0)V n W n . 

Because V* = V- n — > strongly, the second term and its adjoint tend strongly to zero. This 
settles the last assertion. □ 



We can now combine the previous result with Lemma [2.1| and obtain information about 
the strong convergence of the inverses. Note that T(0) is the transpose of T(0). 

Corollary 3.2 Let <fi G L°°(T). IfT n ((f)) is stable, then T((p) is invertible and 

T~\4>) = s-lim T~\4>), T~\4>) = s-lim W n T~\4,)W n . 

If M n ((p) is stable, then M(<p) andT(<p) are invertible and 

M- l ((j)) = s-lim M- x (0), T- l Q>) = s - lim \VnM~ 1 (<j))W n . 
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The next step is the description of invertibility and stability in the case of functions <p 
contained in the set PC of all piecewise continuous functions on the unit circle. For this we 
need more notation. In what follows, let A stand for the set C of continuous functions, the 
set B\, or the set C°° of infinitely diflerentiable functions. We denote by PC [A; K] the set 
of all functions that can be written in the form 

0(e*) = 6(e^)^ + (^)^_(e< e -^)nW( e ' ( ^ r) )V( ei( ^ r) )' ( 4 °) 

r=l 

where b G A is a nonvanishing function with winding number zero, 9±, . . . ,9r G (0, ir) are 
distinct points, the set K C {1,-1, e tdl , . . . , e R , e~ tdl , . . . , e~ t0R } is the set of jump discon- 
tinuities of 0, and (3±, (3f, . . . , (3^ are certain complex parameters. If some these parameters 
are zero, then no jumps occur at the corresponding points. 

The representation (f4Cf ) is essentially the same as (]14|) , however it displays the special 
role of jumps at 1 and —1 and a connection between jumps at a point on the unit circle and 
its complex conjugate. As will be seen below, this distinction is not necessary for the pure 
Toeplitz case, however, it is for the Toeplitz + Hankel case. Note that PC[C; K] is the set 
of all invertible functions in PC with finitely many jumps at K. 

Let PCi[A;K], PCn[A;K] and PCm[A;K] be the sets of functions of the above form 
such that in addition the following conditions (I), (II) and (III), respectively, are satisfied: 

(I) \Ref3±\ < 1/2 and \Ref3±\ < 1/2 for each r; 

(II) -3/4 < Rep+ < 1/4 and -1/4 < Re/5_ < 3/4 and |Re (/?+ + (3~)\ < 1/2 for each r; 

(III) -1/2 < Ref3 + < 1/4 and -1/4 < Re/3„ < 1/2 and 

lRe/3+l < 1/2 and \Ref3~\ < 1/2 and |Re(/?+ + f3;)\ < 1/2 for each r. 

We want to emphasize that PCm[A; K] C PCi[A; K] n PCn[A; K] and that this inclusion 
is proper. This seems strange at first glance, but it finds its solution in the fact that the 
representation (flOP and in particular the /3's need not be unique. Examples of functions 
showing this are given in H, Sect. 4]. Finally, note that <fi G PCi[A;K] if and only if 
4> G PCi[A; K] where 

K = | 1/t : t G K }. (41) 

Invertibility and stability criteria for Toeplitz and Toeplitz + Hankel operators with 
piecewise continuous generating functions can now be stated as follows. For proofs we refer 
to [IIIJ and 0], respectively. 

Proposition 3.3 Let G PC be a function with jumps in a finite set K C T. Then 
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(a) T((j)) is invertible if and only if <p G PCi[C; K}; 

(b) T n ((f>) is stable if and only if (j) G PC\[C] K}; 

(c) M(<p) is invertible if and only if (p G PCn[C; K}; 

(b) M n ((j)) is stable if and only if ' <fi G PCm[C; K]. 

We also need to introduce "approximate" functions </> M , < \i < 1, associated to a 
piecewise continuous function G PC [A; X] of the form (|40|): 

0,(e* e ) = ^ ,e )^,(e ,fl )^,(e ,M )nV.( e!( " r) )V.( e!(e+er) ) ( 42 ) 

r=l 

Here £g ;/1 is the smooth nonvanishing function with winding number zero defined by 

W(e^) = (l - (l - ne~ i9 y P , 0<9<2n. (43) 

Note that also M G A is a nonvanishing function with winding number zero. 

For a (generalized) sequence of operators acting on £2 depending on a parameter 
fi G [0, 1), one can define the concepts of strong convergence (as \i — > 1) and stability in the 
same way as for (discrete) sequences A n . The analogues of Lemma |2j] and Lemma [?72] also 
remain true. We will write A = s — lim A^ if both A^ — > A and A* — > A* strongly. 

Proposition 3.4 Let K be a finite subset ofT, let (ft G PC[C;K], and let M; < \i < 1, 

be the associated approximate functions. Then 

H(cj)) = s-limPT(0 M ), T{4>) = s-limT(^), M(0) = s-limM(^). 

Moreover, if <p G PC[B\;K] and f G C°° vanishes on an open neighborhood of K , then 
H{f4>fji) — > H(fcf)) and H(f/<p /J ) — > H(f /4>) in the trace class norm. 

Proof. First of all note that t p {e ie ) = (1 - e ie ) /3 (l - e -i0 )^, and thus 

^(e^) = exp(2^arg(l-e ie )), W( e *) = exp (2z7?arg(l - /ie ie )), 

where the argument is taken in (— 7r/2, vr/2). It is easy to see — > locally uniformly on 
T \ {1} as fi — > 1. The same holds for the derivatives of arbitrary order. We obtain that 
M — > cf) in measure. Because the sequence M is uniformly bounded in the norm of L°°(T), 
the Laurent operators generated by M (which are unitarily equivalent to multiplication 
operators on L 2 (T)) converge strongly to the Laurent operator generated by 0. This settles 
the first part of the proposition. 
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Now write (p^ = bip^ and (p — bip. From the above statements, it follows that fip^ — > ftp 
in the sense of C°°. Multiplying with b E B\ we obtain that ftp^ — > ftp in the norm of B\ 
and hence the desired convergence of the Hankel operators. The last assertion can be shown 
analogously. □ 

Next we state the necessary and sufficient conditions for the stability of T{(p^) and M(0 M ). 
Proofs are given in jTS and |7[ . These results can be combined with Proposition |3.4j and 
the analogue of Lemma |2J] in order to obtain a result about the strong convergence of the 
inverses and the adjoints of the inverses. We leave the details to the reader. 

Proposition 3.5 Let K be a finite subset ofT, and let (p^, < [l < 1, be the approximate 
functions of the form fc/fy ) associated to a function <p G PC[C; K] of the form ftlQ)- Then 

(a) T(<p[j) is stable if and only if the parameters satisfy condition (I). 

(b) M((pfj) is stable if and only if the parameters satisfy condition (II). 

Note that it does not suffice to require (p G PCi[C;K] or G PCu[C; K], respectively, 
because the representation of <p is not unique and the parameters of M must be chosen 
properly. 

Finally, we need the following basic results. 

Lemma 3.6 Let A : Hi — > H2 and B : H 2 — > H\ be linear bounded operators acting on 
Hilbert spaces H\ and H 2 - Then the following assertions hold. 

(a) The operator I + AB is invertible if and only if so is I + BA. 

(b ) The operator I + AB is a Fredholm operator if and only if so is I + BA. If this is true, 
then ind (J + AB) = ind (J + BA) . 

(c) If A or B is a trace class operator, then det(J + AB) = det(J + BA). 

Proof. Part (a) can be proved by using the formula (/ + BA)~ X = I — B(I + AB)^A. 
Assertion (b) can be proved in the same way, by thinking of the inverses as Fredholm 
regularizers. Also, we use the fact that the kernels (resp. cokernels) of / + AB and / + BA 
have the same dimension. Finally, AB and BA have the same nonzero eigenvalues (taking 
multiplicities into account). □ 

Lemma 3.7 Let A n = P n + P n KP n + W n LW n + C n be a sequence of n x n matrices where 
K and L are trace class operators, and C n tends to zero in the trace class norm. Then 
lim^oo det A n = det(J + K) det(J + L) . 
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Proof. Because W n — > weakly on £ 2 , the sequence KW n L tends to zero in the trace class 
norm. Hence 

detA n = det (/+ P n KP n + W n LW n + C„) 

= det((/ + P n ^P n )(/ + W n LH/ n ) + C;) 

with — » in the trace class norm. Noting that det (J + W n LW n ) = det (I + P n LP n ) 
completes the proof. □ 

4 The limit theorem for piecewise continuous functions 

In Section Q, we have proved that for functions (p G GqB\, the quotient det M n (</>)/ detT n ((p) 
converges to a nonzero constant. Surprisingly, the same turns out to be true for certain 
piecewise continuous functions satisfying a particular condition on the location of the jumps. 
This condition excludes functions with jumps at 1 or —1 or at both a point on the unit circle 
and its complex conjugate. The fact that the Hankel operator H(<f)) need not be compact 
(and hence F(<f>) need not be defined) makes the proof of the limit relation for piecewise 



continuous functions more complicated than the proof of Proposition \Ta . 

In order to overcome this obstacle we introduce another operator determinant which is 
related to F(4>). In this connection, we use the notion of a smooth partition of unity. By 
this we here mean two smooth functions (in C°°) whose sums equals the constant function 
with value one on the unit circle. 

Theorem 4.1 Let G L°°(T). Suppose that T(<f>) is invertible and that there exists a smooth 
partition of unity, f + f = 1, such that H(f<fi) and H(f /</>) are trace class. Moreover, 
introduce the following operators: 

A n = T(/)T- 1 (0)if(0), 

A 12 = T(/)T- 1 (0)P(0)T(7), 

A 21 = T-\<P)H(cf>), 

A 22 = T-\<P)H{ct>)T{f). 

Then An, A 12 and A 22 are trace class, and the operator determinant 

mf) - **(1ZX i + a„ a - A2iAi2 ) m 

is well defined and nonzero. If in addition (f) G B\, then F(<f>; f) = F(<p). 
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Proof. First of all note that if T(0) is invertible, then the function is invertible in L°°(T). 
Using the identities (p3|) and (gj) one can show that An and A 22 are trace class: 



A 22 = T" 1 (0)^(0)T(7) 

= T-\<f>)(H{<l>f)-T{<l>)H{f)^ 

= T-^H^-HiJ), (45) 

An = T(/)T- 1 (0)i/(0) 

= (T(//0)T(0) + H{f/4)H$)T-\<l>)H{i/>) 

= T(f/<f>)H(<f>) + H{fJ<p)H{<p)T~ l ((j))H{(f)) 

= H(f) — H{f /4>)T((j)) + H{f /(j))H{(j))T~ 1 {(j))H{(j)). (46) 

Indeed, in each of these terms there appears a Hankel operator which is trace class. Hence 
also A\2 is trace class. Notice however that A 2 \ need not be trace class. In any case, it 
follows that F(<j>; f) is well defined. Next, we have 



I + A u A l2 \ ( I 0\fI + A n A l2 

-A 21 A U I + A 22 - A 21 A 12 ) { -A 21 I { A 21 1 + A 22 



(47) 



and 



I + A n A l2 \ (I 0\ (T(f) 
A 21 I + A 22 V I J { I 



T~\4>)H{4>){l T{f))- (4f 



Note that F(<fr] /) ^ if and only if the operator appearing in ( |4"7{ ) is invertible, or, equiva- 
lently, if the one in (08|) is invertible. Using Lemma p^§(a) and the fact that 



(/, r ( /)) ( T \f) ) 



it follows that (gSp is invertible if and only if / + T" 1 (0)if(0) = T" 1 (0)M(0) is invertible. 
Hence, we arrive at the conclusion that F((j); /) ^ if and only if M(<j)) is invertible. On 
the other hand, the operator in (f47D equals identity plus a trace class operator, hence it is 
a Fredholm operator with index zero. It follows that so is (|4"S|). Using Lemma |37^(b) we 
obtain that also I + T _1 (0)if(0) and thus M(0) is Fredholm with index zero. Now we need 
only use the fact that M(</>) is invertible if and only if M(<fr) is Fredholm with index zero 
(see [0, Corollary 2.7]). 

Finally, suppose that G B\. Then the operators A u , Ai 2 , A 2 i, A 22 are trace class, and 
one can take the determinant of (|47|) and (fIBD. Noting that the determinant of the first 



matrix on the right hand side of ( |4"T| ) is equal to one and employing Lemma |3.6| (c), it follows 
thatF(0;/)=F(0). □ 
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Although we have not a proof in the general setting, it seems reasonable that F(<p; f) 
does not depend on the particular choice of / (see also the remark below). 

The next result is a quite general version of the limit theorem, which will be applied 
afterwards to certain piecewise continuous functions. 

Theorem 4.2 Let <p £ L°°(T). Suppose that the sequence T n (0) is stable and that there 
exists a smooth partition of unity, f + f — 1, such that H(f<f>), H(f/<p) and H(f/(j>) are 
trace class . Then lim^oo det M n (0)/ detT n (0) = P(0; f). 

Proof. First of all remark that the stability of T n (<f>) implies the invertibility of T(4>). Hence 
the assumptions of Theorem (O] are fulfilled and F(4>] f) is well defined. Let An, . . . , A 2 2 be 



the operators introduced there, and define the following sequences of matrices: 

4? = T n {f)T-\<P)H n {<P\ 

4? = T n {f)T-\<P)H n {<P)T n {f), 

4? = T^(<f>)H n (cf>), 

4? = T- l {4>)H n {<j>)T n {f). 

We first claim that 4i^ ~~ * Ali an d 4s ~~ * ^22 i n the trace class norm as n — > 00. Indeed, 



using fl24|) , fl35|) and (pq), we can write 
4? = T'^HMTM 



T" 1 (0)P„iJ(0)T(/)P n - T?{4>)P n H(<l>)Q n T(f)P n 
T-\<f>)P n (H{<f>f)-T{<f>)H(f))p n -T-\ ( f>)P n V- n H{^^^ 



The last term in the sum tends to zero in the trace class norm because H(<p)H(f) is trace 
class and V- n — > strongly. The first term converges to A22 in the trace class norm since 
T T ^ 1 (0) — ► T _1 (</)) and P* = P n — > I strongly and the expression in the middle is a trace 
class operator. Now we employ identity (|39| ) to rewrite T n (f), and we obtain 

4? = T^T-^HM 

= T n {f/4>)HM + P„P(//0)P(0)P„T- 1 (0)P n (0) 
+ WnHif/^H^WnT-'i^HM. 

Analyzing the first summand yields (see again (P^j), (P3|) and ([?§)) 

T n {f/4>)H n {4>) = P n T{f/<t>)H{4>)P n - P n T{fl4>)Q n H{4>)P n 

= P n (H(f) - P(//0)T(0) N )P n - W n H(f/4>)H{<t>)V n P n . 
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Because H(f /4>)H((f)) is trace class and V* = V^ n — ► strongly, we obtain that 

T n {f/cj>)H n {cj>) - H{f)-H{f/<f>)T$) 

in the trace class norm as n ^ oo. Note that the adjoint of T^ 1 ^) converges strongly to 
the adjoint of T -1 (</>). Hence it is easily seen that 

P n H(f /<P)H(4>)P n T~\<P)H n (<P) - H(f/ (f))H((f))T~ l ((f))H((f)) 

in the trace class norm as n — > oo. Finally, fl55| ) and ([571) imply that 

W n H(f/ ^)H(<p)W n T~ l (<p)H n {(f)) = W n H(fJ^H(<f>)T^W n H n (<i>) 

= H/ nJ ff(//0) J ff(0)T- 1 (0)T(0)K.Pn. 

Because H(f /4>)H((j>) is trace class, the adjoint of T r ^ 1 (0) converges strongly to the adjoint 
of T _1 (0) and V* = V- n — > strongly, it follows that the latter term converges to zero in 

(n) (n) 

the trace class norm. Thus we have proved that A n — > An and A 22 — > ^22 in the trace 
class norm. Now one can immediately conclude that also — > /4 12 in the trace class norm. 
Moreover, it is obvious that A^ — > A 21 strongly. 

The desired limit relation can now be proved as follows. Let 

S n = detM n (0)/detT n (0) = det (P n + T" 1 ^)^)) . 
Since / + / = 1, we have 

p. = (ft, r„(7)) ( T i[> ) . 

This in conjunction with Lemma [JJWc) (for matrices) implies that 



S n = det|p„+(p„, T r Xf))(^ Tn ^^T-\^)H r ^ 

= det I ( P » p) + ( Tn pJ ) T-^HnW (P n , T n {f)) } 



det 



det 



P _L_ A^ n ) 4W 
r n T" /in ^12 

4< n ) p 1 4W 

^21 ■'n ^ ^22 

P 4- 4 (n) 

1 ^11 -^12 

_ An) A n ) p 1 /|( n ) _ /t(™)/t( n ) 
^21 "ll n < ^21 /1 12 
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Note that the last identity follows from 



P \ / p _j_ 4W 4W \ / p 4. 4W 4W 



_/l( n ) p / I ylW p 1 I I /\( n ) A( n ) p 1 /t 00 /I 00 

^21 *n / \ /I21 r n T /io2 / \ ~~ ^21 ^11 "+" ^22 — ^21 ^12 



where the determinant of the first matrix on the left is equal to one. Taking the limit n — > oo, 
we obtain that S n — > /). □ 
Note that if the assumptions of Theorem fO are satisfied (which are slightly stronger 



than those of Theorem |4J]), then the value of F(cp; f) does not depend on the particular 
choice of /. For it is the limit of a sequence independent of /. 

Now we specialize to the case of piecewise continuous functions with jumps on a finite 
set K. We show how the constant F(<ft,f) can be evaluated from the constant F(<f>). The 
proof will reveal that the partition of unity condition is in some sense responsible for the 
afore-mentioned condition on the location of the jumps of 0. Note that this condition can 
be expressed as K n K = where K is defined as in (P|). 

Lemma 4.3 Let <p G PCi[B\;K] with K n K = 0. Moreover, let < p. < 1, be 
the approximate functions associated to cf). Then there exists a function f such that the 
assumptions of Theorem ^TJ and Theorem ^TJ are fulfilled, and F(<p; f) = lim^! F{<p^). 



Proof. The invertibility of T(</>) and stability of T n (4>) follows from Proposition |3.3| (ab). 
Because K C\K — 0, there exists an / G C°° with / + / = 1 such that / vanishes on an open 
neighborhood of K. It is easy to see that ftp G B\ and //</> G B\. Thus H(f<f)), H(f/(j)) 
and H(f/<fi) are trace class. In addition to An, . . . , A 22 defined in Theorem fO] , let 

i00 



AW = r(/)T-(^)F(^), 
A<£ = T{f)T-\<p,)H{<p,)T{f), 



i00 _ 71-1 

L 21 — J 

I 00 _ T"-l 



^21 = T- 1 (0 M ) J E(^), 



^22 = T- 1 (0 M )iJ(0 M )T(/). 



Recall that An and A 2 2 can be written in the form (f45|) and (|46|) , and analogously so can 
and with <\> replaced by 0„. Applying Proposition |3]4] and Proposition |3.5| to these 



modified expressions, we conclude that An —* An and A 2 ^ — > A22 in the trace class norm 
as fj, — > 1. Moreover, also A^ A\ 2 in the trace class norm and A^ — > A 2 % strongly. The 
convergence of these operators implies that F(^>; f) = lim^i Fi^p^; f). Finally, observe that 
^ G G B{. Hence F(<^; /) = F(^) by Theorem 0. □ 
The evaluation of the constant F{<fi\ f) is now straightforward. We consider functions 
<f) G PC\[B\] K] with K fl K = 0. Such a function can be written in the form In terms 
of this representation, the assumptions on the parameters of can be expressed as follows: 
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(i) f3 + = (3- = and \Re/3 r \ < 1/2 for each r; 

(ii) 9 r G (— 7r, 0) U (0, 7r) for each r and 9 r + 9 S for each r and s. 

Hence the final result of this section will nearly treat the situation which was promised in 
the introduction. However, we still must exclude jumps at 1 and —1. 

Corollary 4.4 Let e PCi[B\; K] with K D K = be a function of the form (Qj. Then 

R 

-4B r \ft/2+Pr/2 /, , -w r \P 2 r /2-Pr/2 



F{(f>;f) = F[b]Y[(l-e-^Y r/2+IJr/2 (l + e-^Y 



r=l 

Ft 

r=l l<s<r<B 



Proof. Using Theorem [T5] and Lemma [4.3| we can evaluate F(4>n) and then pass to the limit 
jj, — > 1. It is easy to see that 

log0„(e") = log6( e tf ) + ^A.(log(l-/ic 4 ^)-log(l-A*e^ fl -^); 
Employing the Taylor series expansion of log(l — z) at z = we obtain 

R 

[log0j n = [logfe] n -^^"e-^7n, n>0, 



r=l 



It follows that 



r=l 



= . v ... v R oo 

n=l r=l n=l 

1 oo 

n=l r=l 

i £ /3,/5 s log(l-/. 2 e-^)) 



2^ 1 OT " m 2 

n=l n=l r=l n=l 

oo 

in 

l<r,s<R n=l 
oo R 

-i0 r 

> ii I io'j (;i > i()L» (i_L i //( 

n=l r=l 



2 

Kr,s<i? 



^±(i) \ 1/2 fAilM 1/2 A A N ^ " r/2 



r=l 
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Now it is easy to complete the proof. □ 

Theorem 4.5 Let (f> E PC\[B\; K] with K H K = be a function of the form (TJJ. Then 

detM n (0) ~ G[b] n n nM E M asn-^oo. 

Proof. The previous results say that detM n (0)/detT n (0) — > F(<p; f) where F(<f>;f) is as 
given in Corollary P . Now we need only apply the asymptotic formula for Toeplitz deter- 
minants given in flT5p. □ 



5 The localization theorem for piecewise continuous 
functions 

In this section, we establish a localization theorem for determinants of Toeplitz + Hankel 
matrices. We will show that the quotient det M n (00)/(det M n (0) det M n (if>)) converges to 
some nonzero constant for certain functions and ip. Using this localization, we can reduce 
the asymptotics of detM n (0) for certain piecewise continuous functions to the asymptotics 
for functions of a particular form with only one or two "pure" jumps. 
We first define certain constants in terms of operator determinants. 

Theorem 5.1 Let <j>,xp e L°°(T). 

(a) Suppose thatT(<p) andT(ip) are invertible and H(<f))H(ip) is trace class. Then 

E{<j>^) = det (i + T-^H^H^T- 1 ^ 

= detT- 1 (0)T(#)T- 1 (^) (49) 

is well defined and nonzero. Moreover, E(ip,(p) is well defined and E(ip,(p) = E(<f>,ip). 

(b) Suppose that M(<f>) and M(if>) are invertible and H{ip)M{ip — if)) is trace class. Then 

G{<f),if)) = det (i + M~ 1 ( ( f))H(<p)M(i)-<ip)M- 1 (-il))) 

= det M' 1 ((f)) M^M' 1 (VO (50) 
is well defined and nonzero. 

(c) Suppose thatT((f>) andT(ip) are invertible and H (0) H (if)), H((f))H(if)) and H((f))H(if>) 
are trace class. Then 

H(<f),if>) = E^E&^/E^) (51) 
is well defined and nonzero. Moreover, H(if>, 0) is well defined and H((p, ip) = H(ip, 0). 
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Proof. The equality of the expressions in (|49|) and (|50D follow from (p3|) and (^8|) . Note that 
the operator 

I + T-'i^H^H^T- 1 ^) = T-^^T^T" 1 ^) (52) 

is Fredholm with index zero. Hence so is T((f)ifj). This implies that T(<f>if)) is invertible. It 
follows that also (|52"D is invertible and thus E((f>, ifj) ^ 0. We can argue similarly in the case of 
G(4>, if)) 7^ 0. The point is that the fact that M(cf>ip) is Fredholm with index zero implies the 
invertibility of M(<f>ip) (see [0, Corollary 2.7]). The equality E(<fi,ip) = E(if>,<p) is obtained 
by passing to the transposed operators in (|49|). Recall that the transposed operators of T(0) 
and H(<f)) are T(</>) and H(<f)), respectively. Finally, part (c) follows directly from (a). □ 
Obviously, if G B\, then the above trace class conditions are fulfilled. However, the 
trace class conditions are also fulfilled under weaker assumptions. 

Lemma 5.2 Let (f),ip G L°°(T). 

(a) If there exists a smooth partition of unity, f + g = 1, such that f<j) G B\ and gif) G B\, 
then H(<f>)H(if>) is trace class. 

(b) If there exists a smooth partition of unity, f + g = 1, such that f(f> G B\ and gif) = gif), 
then H(<f>)M(if) — if)) is trace class. 

Proof. Using equation (|24]) we can write 

H{<j>)H$) = H(<P)T{f)H$)+H(<l>)T{g)H$) 

f H{<f>f) - T{<j))H{f))H$) + H{<f>) (H(gi>) - H(g)T(if))) (53) 



Analogously, from equation (p7|), we can conclude 

F(0)M(f-V) = H{4>)T{J)M$-$)+H{4>)T®M$-il>) 
= (H(<l>f)-T(<l>)H(fj)M$-iJ>) 

+H{<f>) (mQ - gifi) - H(g)M(ip - $j) . (54) 



Hence the operators under consideration are trace class. □ 
The next result is a first (general) version of a localization theorem for Toeplitz + Hankel 
determinants. All further work is based on it. 

Theorem 5.3 Let <f>, if) G L°°(T) such that the sequences M n (<f>) and M n (if>) are stable. 
Suppose in addition that H((f))M(if> — if>), H((/))H(if>) and H((f))H(if>) are trace class. Then 

v det M n (<Pifj) ~ ~ 

i^ detM w( 0)detM nW = G ^ E ^ 
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Proof. We obtain from equation ( p8]) that 

and from (|35"D and (|36|) it follows that 

P n M(<f))Q n M^)P n = ^P n T( ( j ) )V n + P n H( ( j ) )V^(v. n T(^)P n + V. n H(^)P r , 

= W n H(0)H(^)W n + W n H(0)H(^)V n P n 

+P n V„ n H(<f>)H(< l p)W n + P n V^ n H{<j>)H^)V n P n . 

Because V* = V_ n — > strongly, it is easy to see that the last three terms tend to zero in 
the trace class norm. Hence 

M n (<^) = M„(0)M„(^) + P n H(<P)M(^ - ^)P n + W n H($)H(if>)W n + C n , 

where C n — > in the trace class norm. Multiplying with the inverses of M n ((p) and M n (if)) 
and applying Proposition jO|(cd) and Corollary |3~2"1 , it follows that 

M-\<P)M n {H)M-\i;) = P n + P n M-\<P)P n H{<P)M{4,-ii)P n M-\i})P n 

W„(^M- 1 (0)W n )^(0)^)(^„M- 1 (^)^„)^„ + C'; 

= J R„ + P n M- 1 (0) J ff(0)M(^-^)M- 1 (^)P„ 
+W n T- 1 (0)H(0)H(^)T- 1 (^)W n + C" n 

with C' n — > and — > in the trace class norm. Lemma |3.7| completes the proof. □ 
Note that the trace class assumptions required in the theorem can be replaced by the 



conditions given in Lemma 5.2 



Now we proceed with establishing basic properties of the operator determinants E(*, *) 
and G(*, *) in the case of smooth functions. These properties allow their computation. Note 
that the constant E(*, *) is already known for a long time ||. 

Theorem 5.4 Let b, c G G Q B\. Then 

E{b,c) = F{b)F{c)/F{bc), (55) 

G(b,c) = E{b,c)/E{b,Z). (56) 

In particular, we have 

E(b,c) = exp ^n[log6]„[logc]_ n V (57) 
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Proof. We start with proving flSSp . Suppose first that 6 = 6 and c = c. Because of (|2D| ) we 
have M(6c) = M(b)M(c). Using the definition ( p3|) of F(*), we can write 

F(b)F(c) = detT- 1 (6)M(6)detT- 1 (c)M(c) = det T _1 (6)M(6) det M(c)T _1 (c) 
= det T- 1 (&)M(&)M(c)T- 1 (c) = det T- 1 (6)M(6c)T- 1 (c) 
= detT- 1 (c)T- 1 (6)T(6c) det T-\bc)M(bc) = E{b,c)F{bc). 

Hence F(b)F(c) = E(b,c)F(bc). Now consider arbitrary b and c. Let b = 6_G[6]6+ and 
c = c_G[c]c+ be the Wiener-Hopf factorization. In the proof of Theorem |2.5| , we have shown 
that F(b) = F(6+6+). Analogously, F(c) = F(c_c_) and F(bc) = F(6+6+c_cL). Doing a 
similar computation as in (p?4T), we obtain 

T-\b)H{b)H{c)T-\c) = T-\b + )H{b + )H{c-)T-\c„). 

Hence E(b,c) = E(b + ,cJ) = £'(6+6+, c_c_). Because 6+6+ and c_c_ are even, we can apply 
the above results. It follows E(b,c) = E(b + b + , cJcJ) = F(6+6+)F(c_cl)/F(6 + 6 + c_c_) = 
F{b)F{c)/ F{bc). Note that ( |55| ) implies fl5T|) by using Theorem |2~5|. 
We are now going to prove (|56|) . We write 

G{b,c) = detM- 1 (6)M(6c)M- 1 (c) 

= det M" 1 (6)T(6) det T- 1 (6)M(6c)T- 1 (c) det T(c)M-\c) 

= det T- 1 (c)T- 1 (6)M(6c)/(F(6)F(c)) 

= det T- 1 (c)T- 1 (6)T(6c) det T- 1 (6c)M(6c)/(F(6)F(c)) 

= E(b,c)F(bc)/(F(b)F(c)) = E(b,c)/E(b,c). 

Here we have used the equations (^), (^), (0) an< ^ (|55|)- 1=1 
Next we address the question under which conditions E(*, *) and G(*, *) are well defined 
for certain piecewise continuous functions and how to evaluate them. 

Lemma 5.5 Let K and L be finite subsets ofT. 

(a) Let <p e PCi[B{\ K] and ip e PCi[B{\ L] with KnL = 0, and let and Vv, < /i < I, 
be the approximate functions associated to (ft and ip. Then E(<fi, ip) is well defined, and 
E((f),ip) = lim hm S(0 w ,V M2 ) = lim lim E{(j)^ip^). 

(b) Let cf) G PCu [B I ; K] and ip £ PCu[B\] L\, assume that there exists an open neigh- 
borhood U of K such that ip\u = ip\u; an d ^ <j>n an d ip^, < /i < \, be the 
approximate functions associated to <fi and ip. Then G((p,ip) is well defined, and 
G(0,<0)= lim lim G(0 Ml ,^ 2 ). 
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Proof, (a) Because K and L are disjoint, there exists a smooth partition of unity, / + g = 1, 
such that / vanishes identically on an open neighborhood of K and g vanishes on an open 



neighborhood of L. With this partition the assumptions of Lemma |5.2| (a) are fulfilled, and 
hence H(<f>)H(if)) is trace class. We write this operator and H^^jH^ip^) in the form (|5*3|), 
and conclude from Proposition |3.4| that if (</> w )if (^.J — > H(<f))H(ifj) in the trace class norm 
as Hi — > 1 and fi 2 —>■ 1. Finally, we apply Proposition |3.5| (a). Note that the order of \i\ and 
/i 2 in the limit does not play a role. 

(b) We choose a smooth partition of unity, / + g = 1, such that / vanishes identically 
on an open neighborhood of K and g vanishes on T \ U. Because f<p G B\ and gip = gip, 
Lemma |5\2|(b) can be applied, and thus the constant G(<fi,ip) is well defined. As if(0 Ml ) is 
trace class, we conclude from Proposition |3.4j and Proposition |3.5| (b) that for fixed /xi 



in the trace class norm as pi 2 — > 1. Writing i7(0 w )M("0 — ip) in the form (|54), we obtain 



from Proposition |3 .4| and Proposition |3.5|(b) that 



in the trace class norm as /xi — > 1. This completes the proof. □ 
The previous result reduces the evaluation of the constants for piecewise continuous 
functions to the case of smooth functions. However, we must warn the reader to be careful 
with the choice of M and ip^. (see the remark after Proposition |3.5| ). It may happen, for 
instance, that given a function G PC\[B\\ K] H PCu{B\]K\, the associated approximate 
functions which satisfy (I) or (II), respectively, are not the same. Later on, we will restrict 
to condition (III), and then this nonuniqueness does not occur. 
Now we establish further relations between these constants. 

Lemma 5.6 Let K, L and M be finite subsets ofT. 

(a) Let G PCi[B\\ K], G PCi[B\\ L], ^ £ PCi[B\\ M], and suppose Kf](LUM) = 
L n M = 0. Then E((f>, ipiip?) = E^ipjE^,^)- If in addition, Kn{LU M) = 0, 
then H{<t>^ 2 ) = H{4)MH{(j)M. 

(b) Let (p G PC m [Bl;K] and ip G PC m [B{;L] with K n L = K nZ = 0, and suppose 
that there exists an open neighborhood U of K such that ip\u = ip\u- Then G(<p,ip) = 

Proof, (a) It follows from (p7|) that the first relation holds for smooth 0, ipi and ip 2 - The 
general case can is obtained from Lemma |5.5| (a) by approximation. Note that ipiip 2 G 
PC\[B\]L U M] and {ip\ip 2 )^ = ipi^fa^ because of L fl M = 0. Finally, the relation for 
H(*, *) follows from the definition by using the relation for E(*, *). 
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(b) The equality for smooth functions is stated in (|56|). In the general case, we ap- 
proximate by smooth functions as indicated in Lemma |5.5| (ab). Note that the approximate 
functions 0^ and M , respectively, are the same for G(*, *) and E(*, *). □ 

Now we establish a first version of the localization theorem for the piecewise continuous 
functions. 

Corollary 5.7 Let cp £ PCui[B\;K\ andip £ PCm\B\; L] with K and L being finite subsets 
of T such that KnL = KnL = ®. Suppose also that there exists an open neighborhood U 
of K such that i/j\u = ip\u- Then 

n-oo det M n (0) det M n (i)) w ^' 



Proof. Similar as in the proof of Lemma |5.5| (see also Lemma |5.2|) one can show that the 



conditions on K, L and U imply that the operators H(<j))H(ifj), H(<ft)H(ifj) and H(<ft)M(tp — 
ijj) are trace class. Moreover, because of Proposition |3.3|, the sequences M n (<f>) and M n (i/)) 



are stable. Hence all the assumptions required in Theorem 5.3 are fulfilled. Note that 



H(<j), i[>) = G(<j>, ifj)E(ijj, 0) by Lemma gj(b). □ 
The previous result requires a "strange" assumption, namely, ip\u = ip\u- Note that in 
case K — (i.e. £ GqB\), one can choose [7 = 0, and hence this assumption is redundant. 
In fact, it turns out that this requirement can be dropped also in general, as shown in the 
following lemma. For technical reasons we sharpen the smoothness condition. 

Lemma 5.8 Let K and L be finite subsets of T with Kr\L = KC\L = $, and suppose 
that ip £ PC[C°°; L). Then we can factor ip = ipiip2 such that ip 1 £ C°°(T) is nonvanishing 
and has ^winding number zero, and that there exists an open neighborhood U of K such that 

■02 1 [7 = ^l\\J- 

Proof. Let U and V be open and disjoint neighborhoods of K and L, respectively, which are 
both non-empty, consist of a finite union of open subarcs, and satisfy U = U and V = V. We 
put ip\{t) = ip(t) for t £ U, and continue ipi on T \ U such that ipi £ C°°(T) is nonvanishing 
and has winding number zero. This is possible because ip is infinitely differentiable and 
nonzero on T \ I D T \ V D clos U. The winding number condition can be fulfilled by 
choosing the values of ip appropriately on some subarc of V. Finally put ip2 — 0/01- Note 
that ip 2 \u = V^lrj = 1- D 

Corollary 5.9 Let £ PCiii[C°°] K] and iji £ PCni[C°°; L] with K and L being finite 
subsets of T satisfying KnL = KnL = ®. Then 

n-oo det M n (<p) det M n (ip) VV ^' ^ J 
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Proof. We factor ip = iftiife as stated in the Lemma |5.8| , and write 

detM n (0^) detM n (0ViV>2) detM n (<^ 2 ) 

X 



detM n (0)detM n (^) det Mj$\) detM n (0^) det M n (0) det M n (^ 2 ) 

detM n (^ 1 )detM n (V' 2 ) 



x- 



detM n (^ 2 



Because of the conditions on K and L and on the function Proposition |3.3| shows that 
the stability of M n (0) and M n (ifj) implies the stability of all other sequences M n (*) which 
occur above. Employing Corollary [5]7| one can take the limit of the above expression, which 
equals 

#(0>l,0-0 2 )#(0,0> 2 ) tt, , i\ rr / / ,x n 

Note that Lemma FTGKa) and Theorem 15. life) has been used here. □ 



Corollary 5.10 Let b G G Q B\, <j)(e ie ) = t Pr (e i{e ~ e ^) and ^(e ie ) = tfs s (e i{9 ' 9a) )- Assume that 
\Rep r \ < 1/2, |Re/3 s | < 1/2, 6 r ,6 s G (-vr,vr], 6 r ^ 6 S and 9 r + 9 S ^ 0. Then 

H(b,tl>) = b + (e ies )H„(e ies )-H + (e~ ies Y% 

H(d>, if>) = (l- e^-^Y rPs (l - e^'A™' (l - e-W^ 



Here b± are the functions defined as in (||) and 

Proof. The calculation is similar to the one given in the proof of Corollary |4.4| . We are using 
(p7|) and Lemma p3|(a). The functions and ip are approximated by M and ip^. We obtain 



E{b,i/>) = b + {e i0s f% 
£(0,00 = (l-e J 



5— Or) 



Because of t p (e- ie ) = t_p(e ie ), we have 0(e ie ) = t^ r (e< e+d ^) and 0>(e* 9 ) = (e^ +e ^). 
The values of the constant H (*, *) follow now immediately. □ 
Now we establish the main result of this section. We consider functions G PCm[B\] K] 
of the form ([II]) and localize as much as possible. This allows us to eliminate functions with 
pure jumps at 1 or at —1. However, using this localization technique one cannot separate 
singularities at both a point on the unit circle and its complex conjugate. 

Theorem 5.11 Let be a function of the form 

R 

0(e*») = b(e w )<J>+(e* e )r(e te )l[Me W ), (58) 

r=l 
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where b e GqB\ and 



(e t0 ) = t p+ (e i0 ), 



M*") = t^(e^%,(e^), 1 < r < R. 

Suppose that 6\, . . . , 9r 6 (0, 7r) are distinct numbers and that 

(a) -1/2 < Re/3 + < 1/4 and -1/4 < Re/3_ < 1/2; 

(b) \Rep+\ < 1/2 and |Re/3"| < 1/2 and |Re(/?+ + (5l)\ < 1/2 /or each 1 < r < R. 



Then 



lim detA W> - = H, 



det Af n (6) det M n ((f) + ) det M n (0~) JJ det M n (0 r ) 

r=l 

w/iere 

if = 6 + (l) 2/3 +6_(l)^+6_(-l) 2/3 -6_(-l)^-2 3 ^- 

R 

J] 6 + ( e ie " )^ ++ ^ 6_ (e* p ) 6+ (e~ i0r )^ + ^ 6_ (e"*' ) ~ ftT 

/3 + (/3++2/3") / . x /3 + (2/3++/3") 



X 

r=l 
R 

X 

r=l 

_ / \ . i I ; f . . — i— ' . > . . I / 

-10,, 



r=l 

/ ... .\PrP7+PtP7+PrPt f . fa , a .\P7Pf+PtP7+PrPf 

X J| h_ e »(«r+«.)j M _ e -»(«r+«.)j 

l<r<s<R 

P+Pt+Pr P7+PrPf f .,„ „ ,\ PtPt+P7P7+PtP7 



X 

l<r<s<_R 



J] (i-e^.))— — — i_ e 



-i{0 r -0 s ) 



Proof. Note that belongs to PC m [B\\ K) with KC{1,-1, e i£ \ . . . , e i9R , e~ ie \ e~ idR } 
being the set of jump discontinuities of 0, and consequently, so does each product which 
involves only some of the factors appearing in (58). We first apply Corollary js]?] and eliminate 



the factor 6. This yields a constant term 

R 

H(b,<i> + )H{b,<f>-)]lH(b,<l> r ) 

r=l 

in the asymptotics (see also Lemma |5.6| (a)). The remaining function (i.e. <p with b being 
dropped) is contained even in PCiii[C°°] K), and we may 
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6 Functions with one jump and the main theorem 



As noted in previous sections, using the localization technique and the limit theorem, we have 
reduced the computation of the asymptotics of the determinants for all piecewise continuous 
functions (satisfying appropriate conditions on the size of the jumps) to those with pure 
jumps at 1 or —1, and to those with two jumps at a point on the unit circle and its complex 
conjugate. 

In this section we compute the asymptotics of the corresponding determinants for the 
functions t l3 (e td ) and t /3 (e^ e_7r )) and thus with these examples the promised asymptotic for- 
mula given in the introduction is proved. 

It is interesting to note that we are able to describe the asymptotic behavior of the 
determinants for the above pure function with arbitrary complex parameters (3. Notice that 
we may exclude the cases (3 G Z as they lead to trivial results. 

In computing these examples, and also in the next section where other interesting ex- 
amples are computed, several Cauchy type determinants arise. The next lemma shows how 
to evaluate several of the products that occur in the Cauchy determinants in terms of the 
Barnes G-function and how to then evaluate the asymptotics that arise from the Barnes 
G-function. It will be used several times in this section and the next. 

Lemma 6.1 (a) For each nonnegative integer n and each z Z we have 

G{l + z-n) _ ( 1 ^ w(w -i) / 2 / sin TT* y G{1 - z + n) 



G(l + z) v ' V * / G(l-z) 

(b) If X! + . . . + x R = yi + . . . + y R and uj := x\ + . . . + x\ - y\ - . . . - y\, then 

R rj\ 4-1" -I- n\ 

~ n^' , as n — > oo. (60) 



T-r Ur^l -T X r Tib) 

11 Gil + y r + n) 

(c) For each nonnegative integers n 1; n 2 and n, the following identities hold: 

J] (k-j) = G(l + n), (61) 

0<i<fc<n-l 

TT (Jfc + i + z) - G(2n-l + z )G (i + i )G(l + l)^ 

11 ^^J^ Z ) - n( „ j_ v s nf „ _ 1 , z\n(„ j_ £ W„-l)(n-2+ 2 ) ' ^ 



0<i<fc<n-l 



11 ( z + k i + k 2) - 7T[Z-Z-^Tr—— v (63) 



G(n + z)G(n - \ + f )G(n + |)2(«- 1 )("- 2 +^) ' 

M 

G(z + n 1 )G(z + n 2 y 

TT + = + " + ni)G(1 7nn (-1) W2( " 2 - 1)/2 ( —Y (64) 



0<fei<ni-l 
0<k 2 <n 2 -l 



o<fei<m-i 

0<fc 2 <n 2 -l 
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Here we assume z ^ {0, —1, —2, . . .} in and ffi^j), and z ^ Z in 

Proof, (a) Using the recurrence relation G(l + z) = T(z)G(z) and the well known formula 
T(z)T(l — z) = ?r/ sin7T2;, we can write 



G(l + z - w) 
G(l + z) 



k=l 



k=l 



71 



. 1 ^n(n-i)/2 f sin^y G(l - 2 + n) 



71 



G(l - z) 



(b) In |L2], Proof of Corollary 3.2] it is shown that 

G(l + z + n) ~ a n b z n n z2 l 2 , as n — > oo, 

where a n and 6 n are sequences of positive numbers not depending on z. 

(c) Noting that G(l) = 1, formula (|61|) can be proved as follows: 



n— 2 n— 1 



j=o fc=i+i 



n-2 
3=0 



G(l+n) 



Formula fl62| ) can be proved in the same sort of way, but requires a little more work. First 
write 



n ( k +j+ 

0<j<fc<n-l 



n— 2 n— 1 



n no-+*+ 



n-2 

n 

3=0 



T(j + n + z) 
T(2j + l + z) 



Now we can write this last product as two products, and then apply the duplication formula 
for the Gamma function. We have 



n-2 



n-2 



Y[T(j + n + z)H 



G(2n -1 + z] 



n-2 

n 



7T 



1/2 



T(2j + l + z) ~ G(n + z) llr(j + i + f)r(i + l + |)2 2 ^' 



3=0 j=0 y J ' v ' 3=0 

and using the basic properties of the Barnes function this is equal to 



G(2n - 1 + z)G{\ + § )G(1 + f )tt^ 
G(n + z)G(n - \ + f )G(n + |)2("-i)(«-2+*) ' 
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Formula ( |B3|) can be shown by writing 

In order to obtain (|64]) , we make an index substitution, and then use (|63|) : 
]~f (z + h- k 2 ) = \ l (z + k 1 + k 2 -n 2 + 1) 

CKfcj^ni— 1 0<fci<»i— 1 

0<fc 2 <"2 _1 0<fe 2 <"2 _1 

G(l + ^ + wi)G(l + ^-w 2 ) 
G(l + 2 + -n 2 )G(l + z)' 

Finally, we apply fl59|). □ 
We now prove the asymptotic formula for the special function tp(e %e ). We first note that 
if A is a matrix of Cauchy type, that is, if {aj}™~Q and {&fc}fc=d are sequences of complex 
numbers such that the following matrix is well defined 



.4 



n-l 



J j,k=0 



then det A = p/q where 



V 



II ( a k-aj)(h-bj), 

0<j<k<n-l 

q = } [ {aj + h 



0<j,k<n~l 

We will use this identity in the next theorem and also for many of the examples that follow. 
Theorem 6.2 Let 4>(e ld ) = t^e 10 ) and assume /3 ^ Z. Then 

det M„(0) ~ n- 3p2/2 -^ 2 (2nf /2 2 3p2/2 G(l/2)- 1 G(l/2 - f3)G(l - f3)G(l + (3). 
Moreover, detM n (0) = if and only if (3 e {1/2,3/2, ... ,n - 1/2}. 
Proof. The Fourier coefficients of = t@ are given by 

sin 7T/9 



LWn ir((3-n)' 
Thus the matrices T n {tp) + H n (tp) have j, entry 



(65) 



sin Ti (3 { 1 1 \ sinvr/3 2/3 - 2j - 1 

+ Tr"! : Z = — ■ 15* a no, , : , To — T. — Zl.i ( 66 ) 



Ti \(3-j + k (3 -1- j -k J 7T f3 2 -(3-2(3j+j+f-k-k 2 
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and except for terms that can be factored out of rows, the corresponding determinant is a 
Cauchy determinant of the above form with a,j = (3 2 — (3 — 2(3 j + j + j 2 and bk = —k — k 2 . 
Our remarks above concerning Cauchy determinants show that 

(o\ n n—1 
71 ' 7= o 



where 



and 



p = 11 (k 2 + k-2(3k-f-j + 2(3j){-k 2 -k+j 2 + j) 

0<j<k<n-l 



q = II (-P + j-k)(l-(3 + j + k). 

0<i,fe<n-l 

We can evaluate the p term by first writing 

V = II (k-j) 2 (k + j + l)(k + j + l-2(3)(-l). 

0<j<k<n-l 

Then we apply ( |61)) and fl62"D to find that p is equal to 

n n - l G{l + n) 2 G(2n)G(l)G(l )G(2n - 2(3)G{1 - (3)G(\ - (3) 



\ ) ^ 



2 2(n-i)(n-i-/9) G ( n + \)G{n)G{n + \)G(n + 1 - 2(3)G(n - (3)G(n + \ - (3) ' 
To evaluate the q term use and (|64|) to see that 



g(l - (3 + 2n)G(l - (3)G(1 -(3 + n)G(l + P + n) (n+1)/2 /shitt/^" 
g G(l-/3 + n) 2 G(l-/3)G(l + /3) 1 j V 77 

We write the product in (|67|) as 



n(2/3- 2J -i) = (-i)"2" r( "Ti^ - 

and then simplify and collect terms to obtain 
detM n (^) = 



(68) 



x 



^2 

G(n + l)G(2n)G(|)G(2n - 2/3)G(l - - /3) 

G{n)G{n + \)G(n + 1 - 2/?)G(n - (3)G{n +\-(3) 



G(l - + n)G(l + P) 
G(l-P + 2n)G(l + P + nY 1 ' 
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In the above expression group together the terms with 2n factors, that is consider 

G(2n)G{2n - 2(3) 
G(l-f3 + 2n) ' 

We wish to apply flBC|) so we multiply and divide by G{2n — (3 — 1) to find that the terms 
involving 2n are asymptotic to 

(2nf 2 - 1 G{2n- (3- 1). 

Before we evaluate the rest asymptotically it is convenient to use the duplication formula 
for the Barnes G-function which reads 

G(z)G(z + \) 2 G{z + 1) = G(l) 2 Tt z 2^ 2z2+3z -^G(2z). (70) 
We let z = n — (3/2 — 1/2 in this formula to obtain 

G(2n - (3 - 1) = 7r - n +P/ 2 + 1 / 2 2 2 ™ 2 - 2 ™/ 3 - 5 ™+/ 32 /2+5/3/2+3 

x G(n - | - §)G(n - § ) 2 G(n - § + |)G(±)- 2 

and then use this substitution in our formula for det M n (tp). 

So at this point we have, gathering all terms, that detM n (tg) is asymptotically 

_/3 2 -l„/3/2-l/2 9 3/3 2 /2+/3/2 ^(|X^ ~ /^(j ~ + P) 

T{n+l- (3)G{1 + n)G(l -(3 + n)G{n - f - \)G{n - §) 2 G(n - | + §) 



x 



G(n)G(n + \)G(n + 1 - 2/3)G(n - /3)G(n + | - /3)G(1 + /? + 



r? 



Write r(n + | — /3) = G(n + §—/?)/ G(n + | — (3) and apply flSDD to the above expression in 
a straightforward way to finally arrive at 

det M n (tp) ~ n- 3 ' 32 / 2 - /3 / 2 (27r) /3 / 2 2 3 ' 32 / 2 G(i)- 1 G(l-/3)G(i-/3)G(l + /5). 

This completes the proof of the asymptotic formula. Finally note that the determinant 
vanishes if and only if the p term or the product fl5S|) vanish. □ 
For functions that have jump discontinuities at the point —1, the analogous results of 
above theorem are contained in the following. 

Theorem 6.3 Let <p(e i9 ) = t p (e^ d -^) and assume (3 Z. Then 

det MM ~ n- 3/32 / 2+/3 / 2 (2vr) /3 / 2 2 3/32 / 2 G(3/2)- 1 G(3/2 - (3)G(1 - (3)G(1 + (3). 
Moreover, det M n (<f>) = if and only if (3 G {3/2, 5/2, . . . , n - 1/2}. 
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Proof. The only effect of the move of the discontinuity is that we need to now evaluate the 
determinant of T n (tp) — H n (tp) since the Fourier coefficients change by a factor of (— l) n . 
In the computation in the previous theorem this replaces the factor (2(3 — 2j — 1) which 
appears in the numerator of the matrices ( |SED with (—2k — 1). A simple of check of the 
computation shows that this only changes the first product in the computation, i.e. the term 
r(n+l/2— /3)/r(l/2— (3) appearing in ( 1571 ) and in the formulas afterwards has to be replaced 
by the factor T(n + l/2)/r(l/2). We leave the details to the reader. □ 
With these two theorems we have completed all the parts of pieces that go together to 
prove formula ([7j). For completeness sake we now state the main result with all the necessary 
restrictions on the j3 parameters. The theorem follows directly from the two theorems of this 
section combined with the localization result of the last section and finally the limit theorem 
of Section £|. 

Theorem 6.4 (Main theorem) Let <p be a function of the form 

r=l 

where b G G Q B\ and where 9i,...,9 R G (— ir, 0) U (0, 7r) are distinct numbers satisfying 
8 r + 6 S ^ for each r and s. Assume also that 

(a) -1/2 < Re(3 + < 1/4 and -1/4 < Re/5_ < 1/2; 

(b) |Re/3 r | < 1/2 for each 1 < r < R. 
Then as n — ► oo, 

detM n (0) ~ G[b} n n nM E Ml 
where the constants G[b], Qm and Em are defined as in M), and 
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7 Other interesting examples 

In this section we evaluate the asymptotics of the determinant det M n (4>) for some classes of 
generating functions with two jumps and for a class of functions that have a singularity of a 
different type. The functions we consider here are special cases where certain assumptions 
on the location and the size of the singularities are supposed. Only one class of functions 
considered here and then only for certain values of the parameters is covered by the previous 
theorems. We begin by considering the following two functions cj)^ 1 '^ and (f)^ 2 ' 13 ^: 

^)(e ie ) = t^ 1/2 (e ie )tp +1/2 (e 1 ^), (71) 

^\e ie ) = tp(e ie )t p (e^). (72) 
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These functions have two jumps at 1 and —1. Note that the functions (j)^ 1 '^ with f3 G Z+l/2 
and (j)^ 2 '^ with (3 G Z are up to a constant equal to the functions (j)(e ld ) 



-And 



hGZ. These 

trivial cases may be excluded without loss of generality. We remark also that (f)^'^ admits 
another representation of a similar form: 



( e *) 



-*S+l/2(e J )</3-l/2(' 



,t(fl-*r)i 



(73) 



In the next proposition we evaluate the Fourier coefficients of c^ 1 '^ and <j)( 2 '^ explicitly. 
Proposition 7.1 Le£ (1,/3) and (2i/3) 6e as above. If (3 £ Z + l/2, taen 

COS 7T/3 



IfP£Z, then 



n(P-n/2) 


sin 7T/3 



i/n is odd 
i/n is even. 



(74) 



«,(2,|8) 



i/n is even 
if n is odd. 



n(J3-n/2) 


Proof. Using the definition ([D of the functions tg, it can be easily verified that 

0&fl(e«) = e%_ 1/2 (e 2 ^), 



(75) 



(2,/3) (e ^ ) 



The Fourier series expansion of 0( 1,/3 ) and (p( 2 '^ can be obtained from that of t$-\ii and t^, 
respectively. Recall that the Fourier coefficients of tp has been given in (|6~5|) . □ 
Because of the special form of the Fourier coefficients of <p( 1 '^ and cj)( 2 'P\ the corresponding 
matrices T n (<f>) + H n (<p) also have a particular structure. In fact, it turns out that they can 
be transformed into matrices of Cauchy form. Thus, using similar computations from the 
previous sections the determinants can be computed. 

Proposition 7.2 Let{a j }^' 1 , {aj}™!^ 1 , {bk}™^ 1 and (b k }^~ 1 such that the block matrix 



A 



A u A12 
A21 A22 



(76) 



is well defined, where 



A u = 


( aj 


+ h) 


-1 


j,k 


0<j 


< 




-i; 


< k < 


nil 


-1, 


A12 = 


( aj 


+ h) 


-1 


j,k 


0< j 


< 


mi 


-i; 


< k < 


m 2 


-1, 


A 2 i = 


(a 3 


+ b k ) 


-1 


j,k 




< 


m 2 


-i; 


< k < 


nil 


-1, 


A22 = 


(a. 


+ h) 


-1 




o< j 


< 


m 2 


-i; 


< k < 


m 2 


- 1. 
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Then det A = p/q, where 



P = ] j (a k - aj)(b k - bj) ] [ (a k - aj)(b k - bj) ] [ (a k - aj)(b k -b 

0<j<k<mi—l 0<j'<fc<m2-l 0<i<mi— 1 

0<fe<m 2 -l 

^ = n + & *) 1 1 + n + II + ^ 

Q<j,k<mi — l 0<j',fc<m2 — 1 



0<i<mi-l 
0<fc<m 2 -l 



0<j<m 2 -l 
0<fc<mi-l 



The matrix A considered in ([F6|) is also of Cauchy form, and therefore the above follows 
from the standard products that arise in the Cauchy determinants. The reason for writing 
A in block form is only for convenience in regard to what follows shortly. 

Now we are able to establish the first main results of this section. In the following theorem 
note that although the values of the parameters (3 — 1/2 and (3+1/2 never fit the requirements 



of the main theorem (Theorem |6.4j ) for any value of (3, the answer agrees with the results 
from that theorem. This can be seen by a straightforward computation taking into account 
the duplication formula for the Barnes G-function ([7(]) with z — 1/2 — (3. Notice, however, 
that if we take the parameters corresponding to representation ([7j^) instead of ([H]), then the 
asymptotic formula given in Theorem FT<| does not coincide with the following result. 



Theorem 7.3 Let (jfi-'M be as defined in ^J\ ) and assume (3 ^ Z+ 1/2. Then 

&etM n (<p^) ~ n- 1 ^ 2 2^ 2 G(l~2(3)G(l/2 + (3)G(3/2 + (3). 

Moreover, det M n ((f)^'^) = if and only if (3 £ {1, 2, 3, . . .} and n > 2(3 + 1. 

Proof. For n fixed, put mi = m 2 = n/2 if n is even, and put m\ = (n + l)/2 and m<i = 
(n — l)/2 if n is odd. Let a = m\ — m^. Denote by Qj k , < j, k < n — 1, the (j, /c)-entry of 
the matrix T„(0( 1,/3 )) + if n (0( 1,/3 ^). We permute the rows and columns of this matrix in such 
a way that we take first the even and then the odd rows and columns. This rearrangement 
results in a matrix B with the same determinant. This matrix is a 2 x 2 block matrix with 
the same structure as ([76|) and with a size determined by mi and m^- 



B 



( 



Q2j,2k 



j,k 



Q2j,2k+1 
£>2j+l,2fc+l 



(77) 



Q2j+l,2k 

Because the even Fourier coefficients of (jp-'ft vanish, it is not hard to see that 

\ 



B 



/ 


i l J 2j+2k+l 


j,k 


'Am 

V2j-2k-\ 


V 


'Am 

Y2j-2k+l 


j,k 


'Am 

S y 2j+2fc+3 



3,k 
],k 



(7i 



J 
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It follows that 



B = -^(' ° I )A, (79) 



7T 



where A is of the form ( ff6|) with a,- = P — j — 1/2, cij = —(3+j + 1/2, bk = —k and b k = k + 1. 
We obtain that 

det£ = f-^^-Y (-l)^det A, (80) 

where det A = p/q with 

p = n n (*-j) a n (i-2/9+*+jxi+*+j) 

0<j'<fc<mi-l 0<j<fc<»ri2-l 0<j<mi— 1 

0<fc<m 2 -l 

. 2 _, 1 , 2 G(l - 2/3 + n)G(l - 2(3) G(l + n) 



G(l + mi) G(l + m 2 ) Q ^ + _ ^ + ^ Q ^ + + ^ 

G(l + m x )G{l + m 2 )G(l -2(3 + n)G{l + n)G{l - 2(3) . 

[HI) 



G(l - 2(3 + m^Gil - 2(3 + m 2 ) 
Q = II G9-1/2-J-*) I] (-/5 + 3/2 + j + fc) 

0<i,fc<mi — 1 0<j,k<rri2 — 1 



X 



0<?<mi-l 0<7<m 2 -l 
0<fc<m 2 -l 0<Kmi-l 



? G(l/2 - (3 + n + a)G(l/2 - (3) G{3/2 - (3 + n - a)G{3/2 - (3) 
' ' 1 G(l/2- p + mi ) 2 ' G{3/2- (3 + m 2 ) 2 



x 



/ i^ ^ 1 / 2 - P + m 1 )G{3/2 + (3 + m 2 ) . w m2 -i )/2 ( costt/3 
1 ' G(l/2- (3 + a)G(3/2 + (3) 1 ' \ n 

G (3/2 -(3 + m 2 )G(l/2 + (3 + m,) r {m _ 1)/2 f cosn(3 



G(3/2 - (3 - a)G(l/2 + (3) v 7 V * 



(_-Qn(n-l)/2 



cos vr/A n (7(1/2 -(3 + n)G{3/2 - (3 + n) 



7T y G(l/2 + (3)G(3/2 + (3) 



G{l/2 + (3 + m 1 )G{3/2 + (3 + m 2 ) 
X G(l/2- (3 + mi )G(3/2- (3 + m 2 )' [ ' 

In order to obtain these results, we have first applied Proposition [77^, then pulled out a 
factor —1 in the first and the third product of the g-term, and finally computed all these 
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products by means of (|6lD, (|63D , and (^4j). Recall that mi + m 2 
observe that 2m x = n + a and 2m 2 = n — a. Note also that a 
whether n is even or odd. Now we can combine (BUD, (KID and 



det 5 = G(l-2/3)G(l/2 + /?)G(3/2 + /3) 



n and ni\ — m 2 = cr, and 
or a = 1 depending on 
), and it follows that 



G(l - 2/3 + n)G(l + n) 



G(l/2 -/3 + n)G(3/2 - /3 + n) 
G{1 + mx)G{l/2 -j8 + mi) G(l + m 2 )G(3/2 - /3 + m 2 ) 



G(l - 2/? + mi)G(l/2 + /3 + mi) G(l -2(3 + m 2 )G(3/2 + /3 + m 2 ) ' 

Here we have used the fact that n(n — l)/2 + m 2 is always even. We apply (pO), and we can 
conclude that the first fraction in the last expression behaves asymptotically as m 3 _1 / 4 ; the 
second fraction as (n/2)^~ 2/3 and the third one as (n/2)~P~ 213 . This yields the desired limit 
behavior of det M n (4>^). It can be read off from the third product in the p-term that the 
determinant vanishes if and only if 2/3 6 {1, . . . , n — 1}. □ 
The results from the next theorem just as in the previous one agree with our main theorem 
(Theorem |6.4[) and this example is partially covered by this theorem. However, we point out 
that the allowed values of (3 + and (3- while of a special form may not satisfy the conditions 
of the main theorem. As before, the duplication formula for the Barnes G-function shows 
the equality of both asymptotic formulas. 

Theorem 7.4 Let cf)^ 2 ^ be as defined in (\7SJ and assume (3^1*. Then 



detM n (0 (2 ' 



~ n 



-W 2 2^ 2 G(l-2p)G(l + p) 2 . 



Moreover, det M n 



(2,0) ^ 



if and only if (3 € {1/2,3/2,5/2,...} and n > 2/3+1. 



Proof. The proof is similar to the one of the previous theorem. We introduce also the numbers 
mi and m 2 and rearrange the rows and columns of the matrix T n ((p^ 2 '^) + Hn^ 2 ' 13 ^) in the 
same way. We obtain a matrix B that can be written in the form (|77|) . However, because 
now the odd Fourier coefficients vanish, formula 



B 



'Am ' 




'Am 








r'2j+2k+2 


j,k 



'Am 




'Am ' 


i l J 2j+2k+2 




V2j-2k 



must be modified as follows: 
\ 

j,k / 



(83) 



It follows that 



B 



sin 7t/3 



7T 



/ 
—I 



A. 



M) 



where A is of the form ([76]) with aj = (3 — j, aj = —(3 + j + 1, bk — k, bk = —k — 1. Hence 

'sin 7r/3 



detB 



71 



(-l) ma det A 



(85) 
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where detA — p/q with 



V 



If -i k -j) 2 II -(^-J) 2 II (l-2P + k + j)(-l-k-j) 



0<j<k<m 1 -l 0<j<k<m 2 -l 0<j<mi-l 

0<fc<m 2 -l 

(-l) n(n " 1)/2 G(l + mi) 2 G(l + m 2 ) 2 

G(l - 2/3 + n)G{l - 2(3) G(l+n) 



x 



G(l -2(3 + mi)G(l - 2/? + m 2 ) G(l + mi)G(l + m 2 ) 

(n _ 1)/2 G(l + m 1 )G(l + m 2 )G(l - 2/? + n)G(l + n)G(l - 2(3) 
[ ' G{l-2(3 + m 1 )G{l-2(3 + m 2 ) ' { } 

n n i-p+j-v 

0<j',fe<mi-l 0<i,fe<m 2 -l 



X 



II (/3 - 1 - J - A:) [I + + 



0<j<mi-l 0<i<m 2 -l 
0<fe<m 2 -l 0<fc<mi-l 



G(l +p + m 1 )G(l + i(mi _ 1)/2 f sinvr/? 



G(1 + /?)G(1-/?) v ' V 7T 

x G(l+/3 + m 2 )G(l-/3 + m 2 ) ^ ^ m2(m2 _ 2)/2 f sinvr/3 



x (-l) mi 



G(l+/3)G(l-/3) ' ' \ 7T 

G(l -/3 + n) 2 G(l -/3) 2 



IH2 



m 2 



^(l-^ + mO^l-^ + m,) 2 



u^/i^/'^VtM G(l+/3 + m 1 )G(l + /3 + m 2 ) 

1 7T J G(l + /3) 2 G(l-/5 + m 1 )G(l-/? + m 2 )" 1 J 



Here we have again employed Proposition [7.2| , pulled out a factor —1 in each of the products 
of the p-term, and finally evaluated the products by fromula floTf). Combining flg5p, (^3) and 
(0) it follows that 

detE = G (l-2/3) G (l + /5 ) 2 .- G(1 - 2/5 + n)G(1 + n) 



G(l-/? + n) 2 

G(l + m 1 )G(l - (3 + mi) G(l + m 2 )G(l - /? + m 2 ) 



G(l - 2/3 + m x )G(l + /3 + mi) G(l - 2/3 + m 2 )G(l + (3 + m 2 ) ' 

By (|60|), the first fraction in this expression behaves as rfi 1 and the second and third fraction 
as (n/2) _2/3 . This yields the limit behavior of det M n ((j)^ 2 '^). Finally, the third product in 
the p-term shows that the determinant vanishes if and only if 2/3 G {1, . . . , n — 1}. □ 
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Next we consider two further classes of functions 0( 3 '#> and c^ 4 " 8 ); 

^\e ie ) = ^ 1/2 (e«)f w (e^)), (88) 
0(^)( e *>) = tpie^/V)^ 6 ^). (89) 

These functions can be directly obtained from (ffi-'P) and <p^ 2 '^ if one rotates them by —n/2 
on the unit circle. The functions 0( 3 >^> and 0^ 4,/3 ^ have two jumps at i and —%. 

The following theorem relates the determinants generated by 0( 3,/3 ) and <^ 4,/3 ) to those 
computed in the previous two theorems. Note that the corresponding asymptotic formulas, 
which can easily be established, cannot be obtained by piecing together the asymptotics for 
two functions with a single jump at % and —%. This fact clearly indicates the limitations of 
the localization idea. 

Theorem 7.5 Let . . . , 0( 4 '^> be the functions defined above. Then 

det M n (0 (3 ' /3) ) = i CT detM n (0 (1,/3) ), (90) 
det M n (0 (4 ' /3) ) = detM n (0 (2 ' /3) ), (91) 

where a = if n is even and a = 1 if n is odd. 

Proof. First note that the Fourier coefficients are related by 4>n'^ = i n (pn'^ ■ We make the 
same rearrangement of the rows and columns as in the proof of Theorem 7\3 and arrive at 
a matrix B for (ffi->w and a matrix B for 0( 3i/3 ) both being of the form (|78|). We have 

S fdiag^)-- 1 \ /diag(^)-- 1 

V °' diag^)^ 1 ; V diagCi 2 *" 1 )^ 1 

Now we take the determinant. Analogously, ^i 4 '^ = V 1 ^^ ■ After the same modification, 
we obtain a matrix I? for <fi( 2 '^ and a matrix i? for <p^ ,,3 \ which are related by 

S Miag^V 1 \ /diag^)^ 1 

V 0' diag(^+ 2 )-VV^^ diag^ 2 )-" 1 

Taking the determinant completes the proof. □ 
Among all the functions we have considered so far in this section, there is only one non- 
trivial function which is even, i.e. which satisfies <p(t) = 4>(l/t), t G T. This is the function 
0( 3 '^ with p = 0. Note that 

^ 3 >°)(e* e ) = <^ 1 1 (92) 

1 -i n/2<9< 3tt/2. 
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However, there are two more interesting examples that can be done which are not piece- 
wise continuous functions, but which are even and have singularities of Fisher-Hartwig type. 

We begin by introducing the function u a , defined by u a (e lB ) = (2 — 2cos6 l ) a . In what 
follows we assume that Re a > —1/2. We first note that if we define the functions 

77 7 (t) = (l-t)\ £s(t) = (1-1/t) 5 , teT, (93) 

where the branches of £ and rj are chosen so that r/ 7 (0) = £s(oo) = 1 for their analytic 
continuations, then u a = £, a Va and tp = C-/3V/3- m the following proposition we list some facts 
already proven in f|T2 , Sect. 3] which will eventually be used to show how the determinants 
generated by u a can be obtained in terms of our previous computations. 

Proposition 7.6 Let D a n be the n x n diagonal matrix defined by 

D a ,n = diag(/if } ,^ a) ,.--^i-i) 



where 



Then 



{a) r(i + a + j) r(i + 7 )r(i + ( 5) 

* j\r{i + a) ' 7,5 r(i + 5 + 7 ) ' 



where the finite Hankel matrix H n (Tg) is defined by 

n-l 



H n (r s ) 



-(i + j)\T(i + sy 
r(2 + i + j + 5) /iJ=0 



The above identities allow us to reduce the computations of the asymptotics for the 
function u a to those for a function that we have already done. 

Theorem 7.7 Let 0(e iS ) = u a {e ie ) with a ^ Z and Re a > -1/2. Then 

Proof. From the above identities we may write 

T n (tp)-H n (t p ) = TZ^DziJ n (^)(T n (^)D^ n ~PH r Xr^)D^ 
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and 



T n (u a ) + H n (u a ) = r a * a D a * n T n (7} a )D2 a) n(T n (£ a )D a * n + aH n {r a )D^ n ). 
Pulling out to the right D~Z n and D~\ and taking the determinant, this yields with (3 
det(T n (u a ) + H n (u a )) det(T-} a D-} n T n ( Va )D 2 ^ n ) det D~] n 



det(T n (t_ a ) - H n (t- a )) " det(r-iD-lT n ( V - a )) detDzl 



Each of these terms can be computed. The determinants of T n (r] a ) and T n (r]_ a ) are equal 
to one since they are triangular matrices. The other matrices on the right hand side are 
diagonal matrices and thus the last term equals 

f} j\T(l + a + j) 

This product can be simplified using the basic recurrence relation of the Barnes G-function 
and an application of (|60|) . The end result is that asymptotically 



det(T n (w Q ) + H n (u a )) _ ^ 2a 2 G(l + a) 



det(T n (t_ a ) - #„(*_«)) G(l - a)G(l + 2a) 

Finally note that the asymptotics of 

det(T n (t_ a ) - H n {t_ a )) = detM n (t_ a (e^)) 



has been computed in Theorem |6.3| . Collecting all terms gives the desired formula. □ 
It is easy to modify this last theorem to find the asymptotic formula for one last example. 
We consider <p(e ie ) = u a (e i(e - n) ) . This is a change in the location of the singularity /zero to 
the point —1. As in the examples of the previous section this change in the location of the 
singularity only requires a small modification in the proof. 

Theorem 7.8 Let (p(e ld ) = u a (e t{e - n) ) with a ^ Z and Re a > -1/2. Then 

Proof. As before the only effect of the move of the discontinuity is that we need to evaluate 
the determinant of T n (u a ) — H n (u a ). This means the only change in the above computation 
is that the term M n (t_ Q ,(e^ £l_7r ^)) is replaced by M n {t_ a (e ld )). Thus using the results of 
Theorem |6.2| the asymptotic formula is established. □ 
Let us make some final remarks concerning the last two theorems. The assumptions a ^ Z 



have been imposed because of corresponding assumptions in Theorem |6.2| and Theorem fx3 
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However, in the last two theorems these assumptions are redundant. To see this, one has to 
elaborate a bit more on the proofs given in the previous section. In fact, one can establish 
an explicit expression for the determinants. In this expression (using analyticity) the term 
G(l — a) appearing in (^) cancels with the term G(l + (3) appearing in 

Note that the condition Re a > —1/2 is exactly the condition for the integrability of the 
function u a . Also this condition can be weakened. One can replace it with the assumption 
2a ^ { — 1, —2, . . .}. Notice that in this case one has to understand distribution with 



well defined Fourier coefficients. For more details we refer to [|12|, [ll] . 
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